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Abstract: We examine the behavior of the spectral function for the operator in QCD 
in the two regimes where it is possible to make analytical progress; weak coupling, and close 
to a second order QCD phase transition. We determine the behavior of the bulk viscosity 
in each regime. We discuss the problem of analytic continuation of the (lattice) Euclidean 
correlation function to determine the spectral function. In each case the spectral function 
has a narrow peak at small frequency; its shape would be challenging to extract accurately 
from lattice data with error bars. 
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1. Introduction 

Ongoing experiments at RHIC are exploring the QCD plasma, achieving initial tempera- 
tures which exceed the QCD phase transition (or crossover) temperature. Since the ions are 
large compared to the intrinsic QCD scale, they create a "macroscopic" sample of plasma, 
which allows the exploration of QCD at these temperatures in its hydrodynamic regime. 
The experimental program indicates that at central rapidities the plasma indeed behaves 
hydrodynamically, developing rather large radial and elliptic flows Q. Viscous hydrody- 
namic analyses Q only describe the observed elliptic flow well if the viscosity is very small, 
rj/s < 0.3. Radial flow is caused by outward pressure accelerating the outer layers of the 
"fireball" radially. It builds up through the whole history of the fireball expansion until the 
hadrons become too dilute to interact. It is therefore sensitive to the whole time history 
of the plasma after the collision. It is sensitive to the equation of state, since this sets the 
relation between pressure and energy density. It should also be sensitive to bulk viscosity. 
The definition of bulk viscosity is a drop in the pressure, relative to the equilibrium value 
at the same energy density, due to expansion: 



with v the flow velocity and £ the bulk viscosity (and with all variables measured in the 
instantaneous local rest frame). Since the bulk viscosity reduces the outward pressure, it 
lowers the amount of radial flow. Of course, under a single set of experimental conditions 
the bulk viscosity may approximately mimic a modification in the equation of state; but 
comparing collisions of different centralities or of different sized nuclei, for which V-i; will 
vary, it should in principle be possible to isolate bulk viscosity from the equation of state. 

At high temperatures where the coupling is weak and the theory is nearly conformal, 
the bulk viscosity is expected to be small [||, ||, ||]. Near the phase transition or crossover, 
however, it may be appreciable. This has encouraged a recent reanalysis of the bulk 
viscosity in QCD. In particular, Kharzeev and Tuchin have recently argued § that the 
bulk viscosity can be determined from lattice results for the equation of state (particularly 
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- 1 - 



the temperature dependence of the pressure P and energy density e) by using an exact 
sum rule and an Ansatz for the functional form of the spectral function of stress-stress 
correlations. And Meyer f7|] has performed a (pure-glue) lattice study of the Euclidean TjZ 
correlation function, with a view towards an analytic continuation of the lattice data to 
determine the real-time spectral function and in particular the bulk viscosity. 

Using Euclidean data to reconstruct real-time correlation functions is in general ill- 
posed without some assumptions about the shape of the real-time correlation function. 
Therefore we feel that it is useful to learn whatever we can about the spectral function for 
in whatever regimes analytic information is available. We have found two such regimes. 
At weak coupling (high temperature) we can compute the spectral function perturbatively. 
And near a second order phase transition (such as probably occurs in realistic QCD at 
some point in the T—p plane), we can make reliable statements about scaling behaviors 
based on universality arguments. The latter case may be of importance for real-world 
QCD; if the critical point is close enough to the temperature axis, heavy ion collisions may 
explore a near-critical crossover, with long correlation lengths and sensitive T dependence 
of thermodynamical variables. 

As discussed above, the bulk viscosity is defined as a deviation of the pressure P = |T? 
from its equilibrium value due to expansion. The operator which generates expansion of 
the system is also ^Tf, and as Kubo showed, one can treat a slowly expanding system by 
coupling this operator to an external source. The bulk viscosity is determined as the linear 
response of the operator |T? to such an external source ||: 



lim — 



dte 



-iuit 



d 3 x 



kTi(x,t), i2?(0,0) 



(1.2) 



Because the energy J d 3 xT 00 is conserved, it is harmless to shift T\ by the energy or any 
multiple of the energy in the above. Two useful choices are k such that (li + kTg) = 
(so the operator we use has vanishing expectation value) and k = — 1, so the correlator is 
replaced with a correlator of T^. 1 

The purpose of this paper is to analyze the spectral function 



p(u) 



dte 



-VJjt 



d 6 x- 



Wx,t),r;(o,o) 



(1.3) 



both in the weak coupling regime and close to the phase transition point, focusing on low 
frequencies. The bulk viscosity is determined by ( = \ lim w _»o p(uj)juj. Section [2] exam- 
ines the perturbative regime, generalizing the bulk viscosity calculation of |5| to nonzero 
frequencies. We show that p/ui has a peak at zero frequency of height 0(a^T 3 ) and area 
f diup/u> ~ o J J 2 T Al . Both the width and area under the peak contradict the Kharzeev- 
Tuchin study ||. Section || studies the spectral function and bulk viscosity near a second 
order phase transition. Drawing on work by Onuki we argue that the bulk viscosity 
shows a power divergence as the critical point is approached; along the crossover line a 
distance t in the T, p plane from the critical point, £ oc t~ zu+a , with v = .630 the scaling 
exponent for the correlation length, a = .110 the critical exponent for the heat capacity, 
and z ~ 3 the dynamical critical exponent. Again, this corresponds to a sharp peak in p/u 



1- We use [—+++] metric conventions. 
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near zero frequency. This result contradicts Ref. [1C]. We end with a discussion section 
in which we examine the implications of our results for the program of determining the 
spectral function via analytic continuation of Euclidean data. 



2. spectral function at weak coupling 

For simplicity we will consider only pure glue QCD here; the behavior of QCD with 
quarks of negligible mass is more complicated but qualitatively the same. The trace of 
the stress tensor is the generator of dilatations, which are a classical symmetry of the the- 
ory. However this symmetry is broken at the quantum level: under dilatations the action 
S = f d^x-^jTr G^ V G^ LU changes because the inverse gauge coupling \/2g 2 changes by 

— /3/g 4 , with (3 = ^j^g 2 ~ <? 4 the beta function for the coupling g 2 . Therefore Tff can be 
replaced with —* -j^TrG 2 the square of the field strength |jlT| . The correlation function 
Eq. ( |1.3j ) is therefore 4/? 2 /g 4 times the correlation of the Lagrangian density with itself. 
The spectral function is related to the Wightman function by a KMS relation: 

G>M = J dte-^J d 3 xi<T^(0,0)T;(x,t)> conn , G>(u>) = 1 _l_ u/T P^)- (2.1) 

It is convenient to shift the operator Tji by a multiple of Tq such that its expectation value 
vanishes, so that the connected correlation function is the same as the full correlation 
function. In our case this means that we need to work not in terms of Tff but in terms 
of O = - T °(r;}/(T ) ~ Tf + 3v 2 T °, with v s the speed of sound. This shift to the 
operator is 0(g 4 ) and will only be important for some terms in what follows. In fact 
the only difference between using O and using Tff is that the correlation function 
will have an extra delta function strictly at zero frequency and of height equal to the 
energy susceptibility (heat capacity) times ((T^)/(T )) 2 . The Wight man correlator for O 
is missing this delta function; its value at uj = is 2T times the bulk viscosity. 

Perturbatively the leading contribution to this correlation function is that of Figure |l|: 




Figure 1: Graph which dominates G 2 -G 2 correlations at weak coupling. The x symbols represent 
the G 2 operator insertions. 

This gives rise to a leading order contribution (defining 2 Q = (u, 0)) of 

G > (Q) = / ^ SHQ-P-R) G> a (P)G>,(R) 

x {jgTP ■ R-P^R v )(g af3 P ■ R-P a R f3 ) , (2.2) 
2 Capital letters are 4- vectors; lower case are their spatial components. 
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where d\ = N 2 — 1 is the dimension of the group. Note that both gauge boson propagators 
are Wightman (cut) propagators, given by 



G> U (P) = [n b (p°)+l] 2^5(P 2 +m 2 ) £ e /Lt (A)e*(A) , (2.3) 

A 

with A the polarization state; e^P^ 1 = 0. The use of massive dispersion relations 3 accounts 
for forward scattering in the plasma; the plasma mass is m 2 = m 2 - > /2 = g 2 N C T 2 /6. Since 
we are working perturbatively, we will treat -C T 2 the energy scale of typical quasipar- 
ticle excitations and we will ignore longitudinal gluons, which have an exponentially small 



spectral weight for p° ~ T [12|. 



There are two contributions to G > (Q), corresponding to the two ways the conditions 
5 4 {Q — P — R)5{P 2 + m 2 )5{R 2 + m 2 ) can be satisfied at vanishing q. Since \p\ = |r|, we have 
p° = ±q°. The contribution where p° = r° = u>/2 gives rise to a continuous contribution, 
corresponding to a cut in the retarded function: 

n> r / oj \ , ,i2 2(3 2 (g) 2d A uj 4 4 4 

Gout = [M T ) + 1J -^4~ ^ 9 oj . (2.4) 

Here 2d\ is the number of color and spin states of gluons. For this contribution, the 
difference between O and is an 0{g 2 ) correction and can be neglected. The parametric 
behavior is simple to understand: g 4 arises because the trace anomaly makes Tji naturally 
0(g 2 ) and we are computing the correlator of two T^'s; and the uj 4 behavior follows on 
dimensional grounds. Note that for uj < T the behavior changes to g 4 uj 2 T 2 because of the 
statistical functions, and is further modified below to ~ gT. 

The other contribution arises when p° = — r°. This requires uj = 0, and therefore 
corresponds to a delta function in the Wightman function (pole in the retarded function). 
If we evaluate the correlation function for operators without subtracting disconnected 
contributions, we find 

G> olc [inc. disconnected] = A J 2 d A (f3 P 2 ) 2 n b (p°)[l+n b (p )] . (2.5) 

This would vanish were it not for dispersion corrections for hard gluons, mentioned above: 
P 2 = —m 2 = — m^/2 ~ g 2 T 2 . Inserting this estimate and taking p ~ T, one finds that the 
delta-function contribution naively scales as g 8 T 5 5(uj), with the g 8 arising as two powers 
of g 2 from the beta functions and two powers of g 2 because of the dispersion relations. 

Since this extra g 4 suppression may come as a surprise for some readers, we will 
pause to discuss its physical origin. Physically, at zero frequency and momentum the 
(x G^Guv operator is probing thermal excitations without disturbing them. A massless 
gauge excitation has equal E 2 and B 2 , and therefore G^G^ = B 2 — E 2 vanishes for an 
undisturbed, propagating gauge boson. It is only because of plasma dispersion corrections, 
which allow E 2 ^ B 2 , that the cancellation is not exact. Therefore the pole contribution 
"wants" to vanish for two reasons; the smallness of the beta function, contributing P 2 /g 4 ~ 
g 4 , and the smallness of dispersion corrections, yielding (B 2 — E 2 ) 2 ~ {g 2 ) 2 ~ g 4 - 



3 We use [—+++] metric convention 
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However there are some complications in evaluating this contribution. First, j3P 2 = 
—[3m 2 is the same order as (Ttf). Therefore the result above is contaminated by discon- 
nected parts, and we should evaluate the correlation function using the operator O defined 
above. This amounts to the substitution —(3m? — > (v 2 — \)p 2 — (3m 2 . The contribution to 
the Wightman correlator is approximately 



^pole 



2 1 l' C 

S ^9 J,, 



n(p)(l + n(p)) 



1 2 

v e 

3 



„2 | m D 

P 19 2 



p 2 dp 
p 



(2.6) 



The next complication is that this integral is infrared singular, indicating that the "pole" 
contribution to the Wightman function is dominated by soft physics. This happens because 
B 2 and E 2 come further from canceling as one considers more infrared, and therefore more 
dispersion-corrected, excitations. The infrared singular, ~ J dp/p 2 , behavior is cut off at 
the scale p ~ gT, where E 2 deviates strongly from p 2 . A complete treatment of this regime 
requires a detailed analysis using Hard Thermal Loop (HTL) effective theory and we 
will not attempt it here. However we can easily see that the linear divergence, cut off at 
the gT scale, reduces by 1 the power of g appearing in the height of the "delta function" 
contribution, such that J^ T dcoG > (u>) ~ g 7 T 5 . Physically, this is because, in the soft 
p ~ gT region, there is no relation between E 2 and B 2 , so the extra g suppression found 
above is absent for such excitations; however they represent only a g 3 fraction of the energy 
density, leading to a delta function contribution ~ g 3 x g 4 ~ g 7 (the g 4 still arising from 
the square of the beta function). 

The last complication is that the "delta function" we just found is not really of zero 
width; interactions broaden it into a sharp peak. Since the bulk viscosity is determined 
by the height of this peak, we need to determine its shape. Interactions mean that the 
delta-function behavior found above receives corrections. The delta function arises from 
integrating over all (thermal) momenta which can run in the loop in Fig. |l]. Very roughly, 
interactions mean that each particle running in the loop contributes not a delta function 
but a Lorentzian of the same area, with width T set by the large angle or large momentum 
change scattering rate for this particle. The width T is momentum dependent and is 
parametrically T ~ g 4 T 3 /p 2 ( see fl3|] for a discussion of the relevant scattering processes). 
Therefore, although soft p ~ gT particles dominate the area of the peak, it is hard p ~ T 
particles which dominate its height, since they have narrower widths. At frequency g 2 T ^ 
u ^> g A T the dominant p is p ~ T(g 4 T /oj) 1 / 2 . The correlation function is parametrically 



g*T*y/T/, 



w<g 4 T 

g 2 T » u > g 4 T 



(2.7) 



For u> ^ g 2 T contributions from the Landau cut and higher order diagrams cannot be 
neglected. We will not attempt to address this region here. 

We can make this estimate more precise by extending the results on bulk viscosity 
(strictly zero w) to treat u ~ g 4 T using the approach of Ref. [14]. The idea is that 
the correlator G > (u>) for uj <C g 2 T is essentially determined by kinetic theory, that is, by 
solving a Boltzmann equation. For a detailed diagrammatic justification for this fact (in the 
context of scalar field theory) see Q. The Boltzmann equation can be solved variationally 



using the tools developed in Ref. [15 for zero frequency and applied to bulk viscosity in 
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||. Ref. [14] showed how to extend these tools from zero to small frequency in the context 
of current-current correlation functions, and there is no problem doing so for stress tensor 
correlation functions too. The computed shape of the correlator G > (u) ~ Tp(uj)/uj for 
u ~ g 4 T is displayed in Fig. |2[ which also compares the shape of the "peak" in the Tff 
correlation function to the peak in the shape of the Tjj — ifTf-k correlation function relevant 
for shear viscosity. Besides the parametrically large difference in the heights of the peaks 
(already discussed above), the figure shows that the peak in the Tij correlator is much 
narrower, more closely resembling a Lorentzian. This is because, unlike the correlator, 
the correlator is not sensitive to soft excitations. 



0.25 



~ o 



Solid: _yT£T|;)(cj)/(a2T 4 ) (for bulk) 



Dotted: _L<T l8 T I8 )(6j)/(TVa*) (for shear) 
2 
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Figure 2: Wightman correlators at small frequency for the operators Tff (relevant for bulk 
viscosity) and — ^-Tkk (relevant for shear viscosity) for N c — 3 pure glue QCD at weak coupling 
(and setting ttid = 1.5T). The "peak" in the bulk correlation function is wider and has much larger 
"shoulders," as discussed in the text. 

Let us compare this weak-coupling behavior to that claimed in a recent analysis by 
Kharzeev and Tuchin || . They derive what they claim to be an exact sum rule showing that 
the frequency integral of the spectral function is related to the temperature dependence of 
the energy density and pressure: 

Eq. (12) of 1: J°° ^-dw = T 5 ^-^- + 16£ vac . (2.8) 

The vacuum energy contribution can be removed by replacing the spectral function with 
its thermal part. They also state that the perturbative contributions should be subtracted, 
though it is not clear to us whether this refers to the righthand, lefthand, or both sides of 
the relation. 

First we analyze the righthand side. At weak coupling and in the absence of quark 
masses it has long been known that P scales as T with a coefficient which can be expanded 



-6- 



in powers of g: 



P(T) = T 4 x (a + Bg 2 \p] (l + t i n ?*\ + . (2.9) 

The coefficients A, B are known but their precise values are not important to this 
discussion. (In fact the coefficients through order g 6 ln(g) are also known [jl?]]). All that 
matters to us here is the functional form and that the scale dependence of g 2 is set by the 
temperature T, which we have shown explicitly in the above by including the 0(g 4 ) term 
with an explicit log which renders the result fj, independent. Using that e = TdP/dT — P = 
3P + 2T 4 B(3, one easily finds that 

< = d - d 7 = \- 2 ^ + °^)-9\ and L^. = 2B(J~g*, (2.10) 

where the coupling g appearing in the expression for (3 should be understood as being 
evaluated at \x ~ T. The dominant T dependence is from this renormalization dependence 
of g: 

Tde-W Td „ 8B ~ , , 
— —, = IB— (3 =—(3 2 ~ g e . (2.11) 
dT T 4 dT g 2 

Therefore a literal interpretation of the righthand side of Eq. ([T^) leads to a parametrically 
0(g e T 4 ) result. 

On the other hand, we can insert our analytical results for the spectral function into the 



lefthand side of Eq. (2.8). The nearly-delta function peak gives a contribution of order g 7 T 4 
(see the discussion after Eq. fl2.6D , and use that p/uj ~ G > /T). The cut contribution, after 
subtracting the vacuum contribution 4 , is 0(g 4 T 4 ). Therefore a literal interpretation of both 
sides of the equation leads to an inequality. And subtracting the rising cut contribution 
does not help, since it changes the LHS from 0{g A ) to 0(g 7 ), not g 6 . Therefore Kharzeev 
and Tuchin's result can only make sense if some subtraction is implied on both sides of the 
equation, in which case it has no utility in the perturbative regime considered here. 



3. spectral function near the 2'nd order transition point 

There is another regime in which it is possible to say something analytical about the T{f 
spectral function. Close to a second order phase transition, low frequency and momentum 
correlation functions are dominated by long wavelength fluctuations which obey universal 
thermodynamical properties which allow for a scaling analysis. In other words, we can 
use universality arguments to determine the functional form of correlation functions, such 
as their parametric dependence on the difference between the temperature and the equi- 
librium temperature. Furthermore, hydrodynamic arguments make it possible to extend 
the universality predictions to dynamical (unequal time) correlation functions, though this 
requires some additional information about what quantities are conserved pl| . 

It is believed that the phase diagram (in the T-fi plane) for realistic QCD with two light 
(but not massless) and one fairly light quark flavor possesses a first order phase transition 

4 After vacuum subtraction, the leading-order cut contribution decays exponentially. However, subleading 
in g 2 contributions do not. According to S. Caron-Huot Jll|, the leading thermal corrections at large lo 
scale as pr ~ g 2 p VSLC (T /uj) 4 ' both for current-current and stress-stress correlation functions, see also JliJ, E3. 
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Figure 3: Expected phase structure in the T — /i plane for realistic QCD. Along the T axis there 
are no phase transitions, but at finite chemical potential there is a l'st order transition line which 
terminates at an Ising universality class critical point. Criticality can give parametric predictions 
for behavior near this point, roughly in the shaded region. At large chemical potential and small 
temperature there is a l'st order transition associated with nuclei and possibly other transitions 
associated with color superconductivity, which are not important here. 



line which terminates in an Ising universality class second order endpoint [22|, as illustrated 
in Fig. U (see however |23| ) . The plasma generated in very high energy heavy ion collisions 
is expected to follow a trajectory close to the T axis (small fi), which probably misses the 
phase transition line [24] but may nevertheless experience a rather sharp crossover with a 
large correlation length. In principle, intermediate energy heavy ion collisions may explore 
the critical point |^5f| . In any case, even if it cannot be directly probed experimentally, 
it is interesting to consider the vicinity of the critical point to see what general lessons it 
teaches us about Tff correlation functions. Also note that the pure-glue theory is expected 
to have a first-order phase transition point [26| but the transition is very weak with a long 
correlation length [27| and so scaling arguments might be suggestive here as well (and two 
color QCD should have a second order transition in the Ising universality class |26||). 

Near the critical point, some linear combinations of the temperature and chemical 
potential map to the temperature and magnetic field variables of the Ising model. Ordinary 
fluids also have a phase transition between liquid and gas phases with a critical point in the 
Ising universality class. There is also a mapping of the Ising variables to the temperature 
and pressure variables of this system, and therefore a mapping between T, n in QCD and 
T, P in the liquid-gas system, as illustrated in Figure f|. 

We need more information to extend universality arguments to unequal time corre- 
lation functions [21|. The long range fluctuations in the order parameter are coupled to 
microscopic degrees of freedom which should lead to diffusive (Langevin) dynamics for the 
order parameter fluctuations. The order parameter should therefore admit a Ginsburg- 
Landau type description, both thermodynamically and dynamically. This leads as usual 
to universality in the behavior of static correlation functions between all critical systems 
with the same dimensionality and underlying symmetries. But at the level of unequal 
time dynamics, the Ginsburg-Landau description must also include any locally conserved 
quantities. In particular, in QCD the temperature and chemical potential are dual to en- 
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ergy density and baryon number density, which are both conserved. Since T and fi map 
to linear combinations of the Ising variables t and h, QCD behaves like an Ising system 
with a locally conserved energy and magnetization. A local, upward fluctuation in the 
magnetization would bias the fluctuations in the order parameter to be positive in that 
neighborhood. Similarly, an upward fluctuation in the (Ising) energy would bias fluctua- 
tions of either sign to be smaller. These local net modifications would persist as long as the 
density of the responsible conserved quantity remained in that neighborhood. Therefore 
there are correlations in the order parameter which persist as long as the conserved quan- 
tities retain their local values. But conserved quantities cannot relax locally; they have 
to move away, which on large length scales occurs diffusively. And the diffusion of these 
conserved quantities is in turn sensitive to the order parameter fluctuations, leading to a 
coupled problem. Fortunately, there is still a notion of universality; systems with second 
order phase transitions in the same (static) universality class, and which have the same set 
of conserved quantities, will show the same dynamical scaling behavior near the transition 
point [§1|. 

For realistic QCD the conserved quantities are the 4-momentum and baryon number. 
The dynamical universal behavior is therefore the same as for the liquid-gas phase transition 
p8|], since this system also has a conserved energy, momentum, and particle number. 




Figure 4: (Color online) Relation between the Ising model temperature and magnetic field 
directions and the directions in QCD and the liquid-gas system. 



Fortunately Onuki has performed a detailed study of the dynamical critical behavior 
of the liquid-gas system near its critical point, focusing on the behavior of bulk viscosity 
|| . Since the dynamical universalities are the same, this can be directly adapted to the 
QCD case 5 . As one approaches the transition point, the correlation length £ exceeds the 
natural microscopic length scale Z m i C (Z m ic = T~ l in a relativistic setting). The dynamics 
of the order parameter fluctuations on a scale Z~* c < k < are dissipative and slowly 
evolving, with the fluctuations changing on a time scale cx k~ z , with z the dynamical 



It may seem confusing that a refativistic theory, QCD, and a nonrefativistic theory, conventionaf fluids, 
shoufd dispiay the same dynamical criticality. But remember that at the hydrodynamical fevef, alf but 
one degree of freedom in either system evolves dissipatively, so propagation of hydrodynamical quantities 
is slow (nonrelativistic) in either case. The exception is sound waves, which propagate at a fixed velocity 
c s < 1 in each theory and which turn out not to be important to critical behavior J2l|] . 
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critical exponent: 

J d 3 xe ik - x (i>(x, r)V(0, 0)} = x^e"^ , r k ~ fc"* , (3.1) 

with ip the order parameter and x(k) its momentum-dependent susceptibility. In 3D Ising 
systems with liquid-gas dynamical universality, z ~ 3 Equilibration is dominated by 
these slow modes. 

In the Ising model, the free energy F is a function of the reduced temperature t = 
(T — T c )/T c and the magnetic field h. Close to the transition point it behaves as 

F(t,h) = F aonsiag (t,h) +t 2 - a F sing (t/\h 1 /^\) . (3.2) 

Here -F S ing is a contribution to the free energy arising from the long range (near)critical 
fluctuations. Note that at small t these fluctuations give almost no contribution to the 
pressure P = —F. The entropy is S = —dF/dt and E = ST + F. Therefore the heat 
capacity C v = dE/dT behaves as C v = C^nonsing + C^singi -0 , displaying a weak divergence 
as t — > which arises from the long range fluctuations in the order parameter. These 
fluctuations provide no pressure but dominate the heat capacity, leading to a speed of 
sound dP/de ~ 0. The same is true in QCD near the critical point except that the Ising 
"temperature" direction corresponds to a linear combination of temperature and chemical 
potential (so the "heat capacity" referred to here is really a linear combination of heat and 
baryon number capacity). 

A small, rapidly applied compression will not promptly change the long-range corre- 
lations of the order parameter; instead it changes the noncritical degrees of freedom, and 
therefore leads to an instantaneous pressure rise AP ~ Ae typical of relativistic degrees of 
freedom. Then, on a time scale r ~ £ z , the long range fluctuations equilibrate, absorbing 
almost all of Ae (since they dominate the heat capacity) and allowing the pressure to relax 
to the equilibrium value. Therefore one might guess that £ ~ £ z . [In the next few para- 
graphs, dimensions in parametric estimates are to be filled in with the appropriate power 
of the intrinsic scale T.] Along the "crossover line" is the map of the h = 0, t > line in 
the Ising system, this behavior is £ ~ t~ zu , with t the distance in the T, [i plane from the 
critical point and v the critical exponent v ~ 0.630 in the Ising system. 

In fact a more detailed analysis || shows that ( ~ ^z-a/u ^ ^-zu+a^ Onuki has 
given a careful derivation of this result; here we give a simple intuitive explanation of 
why it is true. Consider a small, rapid compression. After an intermediate amount of 
time r satisfying T _1 < r < £ z , all modes with k > t 1 / z have equilibrated; those with 
smaller k (longer wavelength) remain out of equilibrium. The heat capacity represented by 
the equilibrated long-range fluctuations is C vk<T i/ z ~ k a / u ~ T a / zlJ _ Therefore a fraction 
T -a/zu Q £ energy remains in the noncritical degrees of freedom so the pressure is elevated 
with respect to equilibrium by AP ~ Aer~ a / zu . This is true at all times r < ~ t~ zu . 
The bulk viscosity is found by integrating this result over all time; it is dominated by 
r ~ £ z , and gives £ ~ ^z-a/u _ p our i er transforming the r dependence gives the spectral 
function (at low frequencies); the Wightman function for frequency u is dominated by 
times r ~ u)^ 1 and is parametrically G > (oj) ~ uj~ 1+a / zu . Reinserting powers of T, we find 

G > {u) ~T 3 (uj/T)- 1+a / zu [u>C z T l ~ z ]\ ( ~T 3 (£T) z - a / u . (3.3) 
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Once again, we find that the spectral function [divided by frequency] possesses a narrow 
peak for uj < T. Note that the total area under the peak is finite and is essentially 
determined by noncritical physics; however the height of the peak diverges as a power of 
the distance to the critical point: ( ~ ^-z/u+a ag one approaches the critical point along 
the crossover line. 



4. Discussion 



In both of the cases where we can gain analytical insight (weak coupling and close to 
the second order phase transition point), the spectral function p{u)/ui has a peak at low 
frequency. This peak arises because there are degrees of freedom which equilibrate very 
slowly; the width of the peak corresponds to the inverse of the relaxation time towards 
equilibrium. At weak coupling, all degrees of freedom exhibit slow relaxation. In this 
case relaxation is slowest for high momentum quasiparticles, which dominate the height 
of the peak, though low momentum fluctuations dominate the area under the peak. Near 
the phase transition point, it is long range fluctuations in the order parameter which 
equilibrate slowly and give rise to the peak in the spectral function. These fluctuations are 
important to bulk viscosity because they dominate the heat capacity, though they are of 
little importance for shear viscosity since they contribute almost nothing to the pressure. 6 
The bulk viscosity is determined by the dynamical critical exponent which determines the 
critical slowing down of these fluctuations. 

It would be too bold to extrapolate from these two examples to claim that G > or p/u) 
always has such a low frequency peak. For instance, at 2T C , where there are no near-critical 
fluctuations but the coupling is strong, it is quite possible that all degrees of freedom 
equilibrate quickly and the spectral function is smooth near zero frequency. However, 
the most interesting temperature range experimentally is temperatures close to T c , both 
because the plasma almost surely explores such temperatures in real heavy ion collisions 
and because it is the only place where we expect the bulk viscosity to be appreciable. 

The presence of a peak in p{ui)/ui at small u is problematic for the reconstruction of 
the bulk viscosity from unequal Euclidean time correlations measured on the lattice. The 
Euclidean time correlation function 7 

G e (t) = J d 3 x^(T^(x,r)T^(0,0)) conn (4.1) 

is related to the spectral function via the integral relation 

f°° du p(u) w cosher- 1/2T)] 

Ge(t)= —^^K(uj,t), K{u,t) = . L , : ' (4.2) 

J_ 00 2tt uj smh{uj/21 J 

In principle, complete knowledge of Ge(t) as an analytic function allows for the recon- 
struction of p{uj)/uj. However in the lattice context one only has numerical data with error 
bars at a finite number of times r, and some procedure (such as the Maximal Entropy 
Method fl29| ) must be used to reconstruct the spectral function. 



6 Note however that it is believed that shear viscosity also shows a weak divergence at the second order 
transition point for systems in the liquid-gas dynamical universality class; roughly r\ ~ £ -0 - 05 . 

7 It is important to compute the connected correlation function on the lattice-or equivalently to subtract 
off the mean value of each operator, which is equivalent at nonzero r to rescaling the T 00 contribution until 
the expectation value vanishes. We thank Derek Teaney for emphasizing this point to us. 
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Figure 5: Kernel K relating the (Minkowski) spectral function to the (Euclidean) unequal 
r correlation function, as a function of (Minkowski) frequency for some representative values of 
Euclidean time r. The behavior at small lu is almost r independent. 

Figure |5| shows the function K(u, r) as a function of u> for some representative values 
of t. This figure illustrates the challenge of learning about a sharp spike in the spectral 
function: all of the curves take the same value near zero frequency. In other words, each 
r gives degenerate information about a peak at small frequency-the area under the peak- 
with almost no sensitivity to the shape of the peak. Therefore, details about the shape of 
such a peak effectively have to be inserted as part of the fitting procedure. Meyer's recent 
results effectively assume that p(lo)/lo is smooth at small to, which we have just shown 
is not a good assumption, at least near a second order transition or a very weak first order 
transition. 

So what is the best way to estimate £ in QCD near the transition (or crossover) 
point? There are two things which we can (in principle) determine reliably from a lattice 
calculation. One is the correlation length £ for the order parameter ((ijjip) for nearly chiral 
2-flavor QCD). The other is the total area under the small momentum peak in the spectral 
function. This can be determined by measuring G~e(t) on the lattice and fitting for the 
spectral function allowing a narrow peak at the origin, with the area under the peak used 
as a parameter in the fit. Call this area Gj = L, <T p{ui)duj/2'KU. [We expect it to be 
almost as large as Ge(t/2).] Then we know parametrically that 

G 

Q = A^T) z - a,v - (4-3) 

What we do not know is the coefficient A. This coefficient requires dynamical informa- 
tion; it represents the constant in the scaling relation between a wave number k and the 
relaxation time of fluctuations at that wave number, Tr^ = A(T/k) z . To get a phenomeno- 
logical estimate of the bulk viscosity where (£T) is large, we have to make some reasonable 
guess for the coefficient A. We advocate A = 1 and A = tt z as two reasonable choices, 
based on the assumptions that the critical regime begins at the scale T and the scale nT 
respectively. Clearly there is a rather large band of uncertainty in the final determined 
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In summary, we have shown that in both of the cases where analytic methods can be 
brought to bear (weak coupling and close to the second order transition point), the spectral 
function needed to determine the bulk viscosity has a narrow peak at low frequency. At 
weak coupling the peak has a height lim^^o p{u)/uj ~ g A T A and area J cLo p(u) / 'oj ~ g 7 T 5 ; 
near the critical point the area is ~ T 5 and the height diverges as T A {^T) z ~ a / u . This 
behavior complicates the reconstruction of the spectral function from Euclidean data. 
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